Unitarity of Maxwell theory on curved spacetimes in the covariant formalism 
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Quantum field theory in curved spacetime may be denned either through a manifestly unitary 
canonical approach or via the manifestly covariant path integral formalism. For gauge theories, these 
two approaches have produced conflicting results, leading to the question of whether the canonical 
approach is covariant, and whether the path integral approach is unitary. We show the unitarity of 
the covariant U(l) Maxwell theory, defined via the Wick rotation of a Euclidean path integral. We 
begin by gauge-fixing the path integral, taking care with zero modes, large gauge transformations, 
and nontrivial bundles. We find an extra geometric factor in the partition function that has been 
overlooked in previous work, coming from the zero mode of the gauge symmetry, which affects the 
entropy and stress-energy tensor. With this extra factor, the covariant calculation agrees with the 
canonical result for ultrastatic manifolds, and in D = 2. Finally, we argue that if there exists a 
unitary (but not necessarily covariant) canonical formulation, then the covariant formulation must 
also be unitary, even if the two approaches disagree. 

PACS numbers: 11.15.-q, 04.62,+v 



I. INTRODUCTION 

The purpose of this article is to resolve some serious 
confusions about how to quantize gauge fields in curved 
spacetimes. There are claims in the literature that the 
canonical and covariant methods of quantization lead to 
different results from each other [J-Q. Analogous issues 
appear for gauge fields in the presence of dielectric media 
and for linearized gravity in curved spacetime Q- 
This is disturbing because it makes it unclear how to 
construct a theory that is simultaneously unitary and 
covariant. 

Since the Maxwell field is free, it is easy enough to pro- 
mote the equations of motion to operator equations. But 
problems arise when one tries to calculate the dependence 
of the partition function Z on the background geometry. 
The geometry dependence is needed for e.g. calculations 
of the stress-energy tensor, or the geometric entropy §. 
To calculate the partition function one needs to gauge- fix 
the action, which can be done covariantly by introducing 
Faddeev-Popov ghosts. Although these ghosts do not in- 
teract with the gauge potential A^, they do contribute 
to the geometry-dependent part of Z. It is therefore im- 
portant to check that they do not introduce any spurious 
unphysical effects. 

We will focus on free Maxwell fields with gauge group 
U(l), on £>-dimensional, spatially compact, connectecu 



and orientable spacetimes. In section [TT] we will discuss 
the properties of the unfixed Maxwell partition function 
on compact Euclidean spacetimes without boundary. In 
section IIIII we will derive the correct form of the covari- 
ant gauge-fixed theory in the continuum, including the 
correct normalization of the partition function Z, which 
is nontrivial. This establishes the equivalence of the fixed 
and unfixed actions on compact manifolds. 

For the most part, our analysis will recover the tradi- 
tional "vector minus two scalars" Faddeev-Popov action. 
But there are subtleties involving zero modes of the gauge 
symmetry and the gauge-fixing terms. These subtleties 
lead to extra terms in the partition function which have 
not been properly taken into account in previous work. 

In section IIVI we will consider two cases in particular 
where the partition function can be calculated by other 
methods: I) the thermal partition function of pure elec- 
trodynamics on an ultrastatic manifold, and 2) the parti- 
tion function on two-dimensional compact oriented man- 
ifolds. We find that the above-mentioned extra terms 
must be included in the partition function in order to 
recover known results. 

Finally, in section [V] we will discuss Maxwell fields on 
nonstatic spacetimes with boundaries. We will show that 
the unfixed covariant theory is unitary after Wick rota- 
tion to globally hyperbolic Lorentzian manifolds. Our 
analysis shows that the covariant formulation of Maxwell 
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1 The assumption of connectedness simplifies notation, but is 



inessential since the partition function of a disconnected man- 
ifold is just the product of partition functions for each of its 
connected components. 
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fields is local and unitary, even if it turns out not to be 
equivalent to the canonical formulation^ 



II. A LOCAL COVARIANT ACTION 

We start with the Euclidean Maxwell partition func- 
tion, prior to gauge-fixing: 



E 



DA 
Vol(G) ' 



-S[F] 



bundles 1 

where S is the standard Maxwell action given by 



S = J d D x^F^ u F^/4, 



(1) 



(2) 



Vol(G) is the volume of the gauge group G consisting 
of all local gauge transformations defined by 8A = Vq 
where a is a scalar parameter which is periodic in the 
range [0, 2w/q), q being the fundamental unit of charge 
associated with the U(l) gauge symmetry When the 
value of a increases by a nonzero multiple of 2ir/q going 
around a closed curve it represents a large gauge trans- 
formation; otherwise it is a small gauge transformation. 

Because the gauge group is U(l), the integral of A 
around around a closed curve is defined only modulo 
2ix/q. As a result there can exist nontrivial bundles cor- 
responding to magnetic instantons. To sum over these 
nontrivial bundles, one simply sums up e~ s for all pos- 
sible harmonic 2-forms satisfying the Dirac quantization 
conditions. 

The path integral measure DA is defined as follows: If 
one writes out modes of A in an orthonormal basis: 



(3) 



then the path integral measure can be written out explic- 
itly as 



DA = ~[[n A dA n , 



(4) 



where fiA is a factor with dimensions of mass required 
to keep Z dimensionless. Similarly, in order to define 
Vol(G), one needs a measure on the space of gauge trans- 
formations: 



Da = rj /4 da r , 



(5) 



where /x^ is a factor with mass-dimension two0 

The partition function (JXJ) has the advantage of be- 
ing manifestly covariant. It is also manifestly local, in 
the sense that e~ s and the measure are products of lo- 
cal functions of the field variables. Unfortunately, it is 
also slightly ill-defined, because the pure-gauge modes of 
A have infinite fluctuations, while Vol(G) is correspond- 
ingly infinite. Formally, these two infinities cancel each 
other out to get a finite answer: the precise finite value of 
oo/oo can be fixed by cutting off each pure gauge mode 
of A at the same value as the corresponding gauge sym- 
metry mode Vq. The zero mode «o corresponding to 
constant gauge transformations needs no cutoff, because 
it is periodic due to the U(l) nature of the gauge sym- 
metry. 

The reader may wonder whether it is really correct 
to include the zero mode cto, since this transformation 
leaves the A field unchanged. The answer is yes. In its 
current form, e~ s and the measure are local functions 
of the fields, being a product of contributions associated 
with each point of space. But if ao were dropped, this 
would no longer be the case. Vol(G) would be a nonlocal 
function of the geometry. In section [V] we will see that 
this locality property is needed to derive the unitarity of 
the theory. Hence, just as Feynman diagrams which are 
invariant under a permutation group Q have to have their 
amplitude multiplied by a symmetry factor 1/|Q|, so also 
a spacetime history whose fields are invariant under a 
global gauge symmetry must have its amplitude divided 
by the volume of the «o moduli space. Integrating using 
the measure (|S|), one finds that the path integral has an 
overall symmetry factor 



(6) 



Vol(ao) 2-Knly/v' 



where V is the volume of spacetime. This factor affects 
the trace of the stress-energy tensor 



T 



SlnZ 



(7) 



as well as the entropy of a thermal manifold with inverse 
temperature /?: 



(8) 



In section HVl we will show that this factor is necessary in 
order to obtain correct results for thermodynamic calcu- 
lations. 



2 We originally attempted to use BRST to prove unitarity, but it 
turns out that BRST is problematic for Maxwell fields on spa- 
tially compact manifolds. The representation of BRST is re- 
ducible [jj, because there exists a canonically conjugate pair of 
spatially constant ghosts which are unpaired with any unphysi- 
cal vector modes 1101) . One can try to deal with this by adding 
ghosts-for-ghosts, but we were unable to deal with the resulting 
IR divergences in a satisfactory way. 



3 These factors of /i are similar to a UV cutoff A, insofar as they 
are dimensionful parameters needed to make the path integral 
well-defined. Indeed, in some regulator schemes (such as the 
heat kernel), the same dimensionful parameter plays the role of 
both A and fi. Nevertheless, it is important to recognize that 
this (j, does not result from UV divergences, and can therefore 
play a role even for D = 2 Maxwell fields, where there are no 
local degrees of freedom. 
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III. GAUGE-FIXING 

We now perform the Feynman-Faddeev-Popov-t'Hooft 
trick in order to gauge-fix the small gauge transforma- 
tions. We would like to insert into the path integral a 
factor of unity which includes the delta function 



£(V^-cj), 



(9) 



with uj being an auxilliary scalar parameter. However 
this is overconstraining, since on a compact surface one 
of the constraints is redundant (if the integral of uj is 
zero) or impossible (if the integral of uj is nonzero) . 

To deal with this in a covariant way, we use a complete 
basis of orthonormal scalar modes 4> n of the Laplacian: 



omit the large gauge transformations, which will be han- 
dled separately. 

The determinant can be viewed as the quantum field 
theory of two scalar Grassmannian fields known as 
Faddeev-Popov ghosts, whose partition function is given 
by 



det't-^V 2 ) = / D'cD'ce 



where the Grassmannian integrals are defined by 



(13) 



- V 2 </>„ = X n 4>n, 



d D x\/g<t>i, 



l. 



(10) 



This allows us to expand out uj and V m j4 m in orthonormal 
modes of the Laplacian: 

U = J2"n<l>n, V M i4" = 53(V/*i4")„tfn. ( U ) 

n n 

We can then insert into the path integral a factor of unity 
which includes every component of the delta function ex- 
cept the zero mode: 

1 = / D'aU *((V„Af a) )„ - u, n )det'(- M - 2 V 2 ), (12) 

where A^ a) = + V^a, det'(-/i~ 2 V 2 ) is the usual 
Faddeev-Popov determinant, and the primes arc re- 
minders to omit zero modes. In the case of D'a we also 



J 



D'cD'c= 



d d 
dc n dc„ 



(14) 



with the zero modes of c and c omitted. The factor of 
fi in Eq. (|14[) is necessary to make the partition function 
dimensionlcss; the fact that it is the same fx a that appears 
in Eq. ([5]) is required so that the Faddeev-Popov trick is 
equivalent to inserting an identity operator. 

One can then perform the Da integrals by changing 
variables from A to Ai a -\, making the integrand inde- 
pendent of a. leading to a factor of Vol(G) that can- 
cels (formally) with the one in Eq. ([1]). This eliminates 
all gauge symmetry except for the zero mode and large 
gauge transformations, giving the following gauge-fixed 
partition function: 



£ J ^-D'cD'c J] <5((V^)„ - uj n ) e 

bundles large n^O 



-S[F,t 



(15) 



The measure -D^4/Gi argc is an integration over equivalence classes of A under large gauge transformations; formally 
it can be thought of as dividing by the number of Gribov copies, i.e. configurations of A which satisfy V^A^ = and 
are equivalent to A = under a large U(l) gauge-transformation. 
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We can then give dynamics to the unphysical modes 
of A by integrating over uj using the identity 



partition function: 



N(0 



Dwe^' iDl * 2 ^ = 1, 



(16) 



= J T\lrv £ f ^-D'cD'ce- 3 ^, (17) 
M 2 V 2nV ^±Z_J Gi arge 



bundles 1 



where the action is the usual covariant action in t'Hooft 
Re gauge 



where N is an infinite normalization factor. Since this 
normalization factor is local, and depends only on the 
spacetime volume yfg, it can be dropped by absorbing it 
into a redefinition of the cosmological constant. One can 
then perform the Duj integrals to obtain the following 



S 



d D x y /g±F lu ,F> u ' + hv ll A> l )*-cV 2 c, (18) 



and there is an extra factor coming from the Gaussian 
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integral over the unmatched zero mode ujq: 
duj e -"o/ 2 « = 



(19) 



In this gauge-fixed form of the partition function, covari- 
ance is manifest, but locality has been hidden because 
the zero modes of the field are treated on a different foot- 
ing than the nonzero modes. Nevertheless, the theory is 
equivalent by construction to the original unfixed action 

©• 

To further evaluate the partition function of Eq. (fTT|) , 
we divide the fields into the zero modes and non-zero 
modes. The non-zero modes can be evaluated directly as 
Gaussian integrals, 



Zg 



TT ft* 
11 ^VT 



n ^v? 

Aeo-s 

n ^ 2x 



(20) 

(21) 
(22) 



Here Z? is the path integral over transverse modes, which 
depends on cry, the nonzero spectrum of the Hodge 
Laplacian on 1-forms restricted to act on transverse 
modes. Z^ is the path integral over longitudinal modes, 
and Zg is the path integral over ghosts; both of these 
depend on erg, the nonzero spectrum of the scalar Lapla- 
cian. 

Now consider the zero modes of A. The path inte- 
gral over zero modes does not have the form of a Gaus- 
sian integral, as was the case for the non-zero modes. 
The path integral over these modes would be infinite, 
were it not cut off by identifying configurations related 
by large gauge transformations. Large gauge transfor- 
mations identify A and A + (2-Kjq)r\ for r\ G H 1 (M, Z), 
where the cohomology group H 1 (M, Z) is the set of har- 
monic 1-forms whose holonomies are all integer- valued. 
The space of flat connections modulo large gauge trans- 
formations is the moduli space and the path integral 
over zero modes is equal to its volume Vol(^#). 

It will be convenient to adopt a non-orthonormal basis 
for ^# consisting of harmonic 1-form fields w% , . . . , Wb 1 
that form a Z-basis for H 1 (M, Z) , i.e. such that every 
1-form with integer-valued holonomies can be written as 
a linear combination of the wi with integer coefficients. 
In this basis, the flat connections modulo gauge trans- 
formations is the region [0, 2ir/q) bl . The price to pay for 
this choice of basis is that the path integral measure ac- 
quires a factor \J det(L), where Y is the metric on moduli 
space: 



L/,/ = / wi A *wj = J d Xy/gg a wi a Wjb 
The volume of the moduli space is therefore 



Vol(^) 



Vdet(r). 



(23) 



(24) 



Finally, we have the monopole configurations. These 
are indexed by harmonic two-form fields F satisfying the 
Dirac quantization condition, F 6 =^H 2 [M, Z), where 

the cohomology group H 2 (M, Z) is the set of 2-forms 
having integer flux through every closed two-dimensional 
surface ScE. The contribution of these field configura- 
tion to the partition function is 



Z 



bundles 



E 



-S[F] 



(25) 



Thus combining all these factors we obtain 



Z = ^\l7T77Vo\(JZ)Z T Z L Z< 
//' V 2irV 



G 



E ■ 

^# 2 (M,Z) 



-S[F] 



(26) 



IV. EXAMPLES 

A. Two-dimensional case 

Let us consider the case where the manifold M is 
two-dimensional and orientable. In two dimensions, the 
Maxwell field has no local degrees of freedom. This is 
reflected in the fact that the path integral simplifies con- 
siderably, due to a cancellation between the gauge po- 
tential and ghosts. In particular, the nonzero spectrum 
ctt of the transverse Laplacian is identical to the nonzero 
spectrum of the scalar Laplacian 175 , so the contribution 
of the nonzero modes is 



z T z L z G = n 4-^\/£ 



pi 



1 



AG<ts 



(27) 



Here we used the fact that if the product had been a 
product over all scalar modes (not just non-zero modes) 
then it could be absorbed into a local term@ Note that 
this factor of y/£ cancels the one in Eq. (flT)) , as required 
by gauge-invariance of the partition function. 

To compute the volume of the moduli space, we appeal 
to Poincare duality to show that the metric F ([2"3"| on 
moduli space has unit determinant In two dimensions, 
Poincare duality implies that under the bilinear pairing 



<»?,c> 



r] AC, 



(28) 



the lattice ^(M^Z) is self-dual: 

{i):(i|,()e2 VCe H 1 (M,Z)} = H 1 (M,Z). (29) 



4 In this step we are essentially just reintroducing the factor of 
JV(£) that was dropped from Eq. <16t . 

5 For a concise summary of the relevant aspects of Poincare duality 
see Ref. Ill 
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Expressed in the basis wi, this implies that the matrix 

Pu = JwjAwj (30) 

is invertible over the integers, and therefore dct(P) = ±1. 
The metric Yu is related to the bilinear form (|28|) by 

rij = (wi,*wj) (31) 

where * denotes the Hodge star. In terms of matrices, 
r = PS where S is the matrix *wi = Sijwj. Since * 2 = 
— 1, the matrix S satisfies det(S) = ±1. The determinant 
of r is then dct V = dct P dct 5 = 1 (the positive sign 
follows from the fact that Y is positive-definite). Thus in 
two dimensions, the volume of the moduli space is 

voi(^) = Vdclr = (^p)" 1 . (32) 

On a two-dimensional connected manifold, the only 
harmonic two-forms are proportional to the volume form 
e a t>, and the Dirac quantization condition gives 

tp 2n n fi^ 

Fab — —yt-ab [o6) 

where ngZ. Their contribution to the partition function 
is 

^bundles = Y j e- 2 ^ n2 ' V " 2 . (34) 

Combining these factors we obtain 

Since we assume M is connected and orientable, the ex- 
ponent bi — 2 = —x, the Euler characteristic of M. Us- 
ing the Gauss-Bonnet theorem X = J y/ljRi this term 
can be absorbed into a local counterterm (the Einstcin- 
Hilbert term) and we therefore discard it. By Poisson 



summation the remaining part of the partition function 
can be put in the canonical form 

Z = Y J ^ hVq2n2 - (36) 

n 

This agrees with the known result for the partition func- 
tion of two-dimensional gauge theory, calculated using 
topological quantum field theory methods [l2j . 



B. Ultrastatic manifold 

We now consider the ultrastatic case in which M is of 
the form S 1 x E, where the S 1 direction is the imagi- 
nary time coordinate r with period j3 and E is compact 
and orientable, but otherwise arbitrary. We will show 
that the formula ([26]) reproduces the canonical form of 
the partition function Z = tre~@ H . We divide the field 
modes as follows: 

1. Nonharmonic modes of E: 

(a) transverse spatial polarizations; 

(b) longitudinal spatial polarizations; 

(c) temporal polarizations; 

(d) ghosts; 

2. Harmonic modes of E: 

(e) vector zero modes of E; 

(f) spatially constant temporal modes; 

(g) spatially constant ghosts; 

3. Nontrivial U(l) bundles on SxS". 

In the first category, there is a factor of 
(2 sinh(/3-\/A/2)) _1 , the canonical partition function 
of a mode with frequency vA, for each transverse vector 
mode with Laplacian eigenvalue A. The contributions of 
the remaining modes to Z cancel each other exactly. 



The second category includes, not only the spatially constant temporal mode and ghosts, but also b\ — 1 vector 
zero modes of E. Here b\ is the Bctti number of M (which counts the number of harmonic vectors); we subtract 1 
from this because the temporal mode is included separately. Excluding the integral over vectors that are constant in 
time (the moduli space), the contribution of all these modes to Z is 

n 9 / 27m \ 2 r- — - / 27m \ 
*^(-j 

where we have used ^-function regularization of the product, n„>i an = a^^e - ^ (°) = ^2ir/a. Note that the factor 
£ -1 / 2 cancels the explicit factor of in Eq. (f26|) . as required by gauge invariance. 

I 

Next we calculate the volume of the moduli space. A form wi = /3 dr wrapping the time direction together 
basis for the cohomology group of M is given by a one- 



/'27m 
Mav2tt [ — 



(37) 
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with a basis W2, ■ ■ ■ ,Wb 1 of i/ 1 (E,Z). Letting 7 denote 
the metric on the moduli space of E, the metric T on the 
moduli space of M is given by 



r 



and hence its volume is 







hu 



(38) 



which is precisely the contribution of the magnetic 
monopolcs to the canonical partition function. 

Now consider the electric fields: 2-forms proportional 
to drAwi . If we integrate over a surface S 1 x 7/ the quan- 
tization condition implies that F = ^njdr A wj where 
nj G Z. The contribution of these field configurations to 
the partition function is 



Vol(^T) = (^p) bl /3 (bl ~ 2)/ V^ det 7 . (39) 

The final kind of contribution comes from nontrivial 
U(l) bundles. These are proportional to suitably quan- 
tized harmonic 2-forms on M . They can be divided into 
two types: harmonic 2-forms on E constant in r (mag- 
netic fields) , and those of the form dr A wi where wi is 
a harmonic 1-form on E (electric fields). 

The magnetic fields are spatial 2-forms B satisfying 
the Dirac quantization condition J s B £ 2^Z for every 
closed two-dimensional surface S. These contribute to 
the partition function as 



2 jj 



(40) 



E 



exp 



--(-)y 

2/3 \q I Z - 



Jijninj 



(41) 



This is not of the canonical form, because the factor of 
f3 in the exponent of Eq. (|4ip appears in the denomi- 
nator. However, we can remedy this using the Poisson 
summation formula, 



-1/2 



e~ xl ' Mx = det ( — ) e"^ TM ~\ (42) 



where M = jp(^-) 2 J- Thus we can rewrite Eq. (|4*Tj) as 



Z E = 



2ir 



(6i-l)/2 



V / det(7) 



£ exp 



59 2 /j 
—7 mimj 



(43) 



where 7 is the inverse of jjj and Einstein summation is assumed. 



The prefactors in Eqs. (126J) , ((37J) , (J39J) , (J43J) cancel, so the 
path integral has the form of a canonical partition func- 
tion. All that remains to show is that the sum in Eq. (|4"3")l 
is the same as the canonical partition function of the 
harmonic i?-fields satisfying the quantization condition 
J s E E ql for all closed (D — 2)-surfaces S. In what fol- 
lows we treat E as an (D — 2)-form related to the usual 
vector field E by Hodge duality of vector fields on E. By 
Poincare duality, the condition that J s E € qL for closed 
S is equivalent to J wi AE = qmj where mj € Z for all I. 
Since E is a harmonic (D — 2)-form, it can be expanded 
in the basis {*u>/}, and takes the form 



E 



q^y IJ mi * Wj. 



The energy of such an E'-field is 

1 

2 



EA-kE 



n 2 

1 IJ 

y7 m/TOj. 



(44) 



(45) 



Thus the sum in Eq (|43p is precisely the canonical par- 
tition function of the harmonic E'-fields satisfying the 
quantization condition. 

The path integral expression (|26|) therefore gives the 
same result as the canonical partition function. Hence 
the covariant formulation of quantum electrodynamics is 
unitary on ultrastatic manifolds. 



I 

Note that this result depends crucially on using the 
correct normalization factor in the gauge-fixed partition 
function [251 Without this prefactor, we would not get 
the right result even in the trivial case of D = 1 electro- 
magnetism, where (because E is 0-dimensional) there are 
no physical degrees of freedom at all. Yet the Faddeev- 
Popov action has two ghosts but only one vector mode, 
so it appears at first sight that there are -1 scalar degrees 
of freedom. This spurious contribution is cancelled when 
the zero modes are properly handled. 



V. NON-STATIC UNITARITY 

On non-static manifolds, it is not clear that the co- 
variant and canonical approaches agree with one another 
Q. Nevertheless, we will argue that the covariant parti- 
tion function must be unitary once it is Wick rotated to 
Lorentzian signature (assuming without proof that this 
Wick rotation is possible) . Our argument will rely on the 
existence of a unitary canonical approach to clectromag- 
nctism, but we will not assume that the two approaches 
agree in all respects. 

In this section, we will go back to the original unfixed 
form |T]) of the covariant Euclidean partition function Z. 
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In order to define unitarity, we need to be able to evaluate 
Z on a manifold M with boundary dM . 

A. Manifolds with boundary 

States. The partition function of a manifold M with 
boundary dM is a functional of the geometry of M and 
of a specified boundary value of the dynamical holds. Be- 
cause the unhxed action does not propagate the temporal 
or longitudinal modes, in order to compute Z(M) it is 
sufficient to specify the connection on dM. Let Aqm re- 
fer to the pullback of the vector potential to dM (which 
includes only the D — 1 components of A which point 
along dM). Then the connection is specified by Aqm 
modulo gauge transformations cxqm- 

Since the value of Z(M) depends on Ag M , the path 
integral can be regarded as outputting a wavefunction 
^[^4-3m] on the space of boundary vector potentials. This 
wavefunction can be written as 

bugles J V01 (Ga/) 

where the integral is over connections with a specified 
pullback to dM and Gm refers to the group of gauge 
transformations a which vanish on dM. However, gauge 
symmetry implies that is constant along modes of A 
which arc pure gauge on dM. 

Inner Product. For any region E C dM, the inner 
product on the space of states is given by 

<*,*>== E / vorfe) $ * [j4s] *^ E] ' (47) 

bundles 

where Gs refers to gauge transformations on E that van- 
ish on <9E. Any two open manifolds Mi and M 2 can be 
glued together along a shared boundary E C dM\ , dM 2 
by using the symmetric bilinear (<3?*, to evaluate the 
pair of states at E. In other words, one requires the fields 
at the boundary to match, and integrates over all pos- 
sible fields while being careful to mod out by the newly 
introduced gauge symmetries. By doing this, one obtains 
Z for the combined region Mi Us M 2 , where Lix denotes 
the union of two disjoint submanifolds, glued together at 
a shared boundary A0 

Superselection Sectors. Suppose that a spatial slice 
E is decomposed into two open regions joined by a bound- 
ary B, so that E = Ei Us E 2 . We would like to be able 
to say that the Hilbert space decomposes into ten- 
sor factors T-Ls 1 ® "Hs 2 . However, this is not quite true 



This subscript will sometimes be dropped when the identity of X 
is unimportant to the discussion, but the step where the mani- 
folds arc glued together at their shared boundary is still implied. 



since the normal components of the electric and mag- 
netic fields must match across the boundary B. Since 
these fields can be measured on both Ei and E 2 , by mi- 
crocausality they must commute with all other operators 
localized in either Ei or E 2 , i.e. from the perspective of 
a single region they are superselection sectors. 

In order to implement the magnetic constraint, we re- 
quire that 1) the choice of bundle on B must match on 
both sides, and 2) so must the connection; in other words 
the vector potential Ab must agree up to gauge trans- 
formations (including large ones). Note that if B has 
noncontractable curves, this constraint also requires the 
holonomy of A around those curves to agree. 

The electric constraint has to do with Gb, the group 
of gauge transformations on B. The fields in a single 
region (e.g. Ei) do not have to be invariant under gauge 
transformations which affect the boundary. Instead they 
can transform under nontrivial unitary representations of 
Gb- Since Gb is abelian, the irreps are all 1-dimcnsional, 
and are labelled by the choice of electric flux across B. 
However, the state of E must be invariant under Gb, so 
the irrep of Ei must be conjugate to the irrep of E 2 . This 
enforces the constraint that the normal component of the 
electric field matches across the boundary. 

Evolution Operator. Now suppose we divide the 
boundary of a Euclidean spacetime region into two parts 
joined by B, so that dM = Ei U_g E 2 . Then the partition 
function can also be used to calculate a time evolution 
operator £ cov (Ei — > E 2 ) which evolves a state of the fields 
in Ei, to a state of the fields in the complement E 2 , while 
preserving superselection data at B: 

£ cav = (gJ®I)Z(M), (48) 

where J is the antiunitary time-reversal operator (i.e. 
complex conjugation of the position-space wavefunction), 
and g is the inner product defined in Eq. (|47[) . Here 
Z(M) is regarded as a state using Eq. (|46|) . gj is an op- 
erator whose domain is the Hilbert space but whose 
target is the dual vector space / H^ 1 , and / is the identity 
acting on Hs 2 . Thus £ cov takes states from "Hsi to Hs 2 . 

Wick Rotation. We now wish to define £ cov for a 
globally hyperbolic Lorentzian manifold M with bound- 
ary dM = Ei Ub E 2 , such that Ei is an initial spacelike 
slice and E 2 is a final spacelike slice. To do this, we con- 
sider the class of all Euclidean metrics for which dM has 
the same geometry. £ cov is defined in terms of the parti- 
tion function Z(M), which in turn depends on the metric 
9nw We will assume without proof that, as usual in QFT, 
Z(M) is an analytic function of g^, which can be ana- 
lytically continued through complex values of g^ v to the 
Lorentzian signature, without encountering obstructions 
or ambiguities. Because the Euclidean partition function 
is real, the Lorentzian partition function is invariant un- 
der time-reversal T. 

Converting from Z(M) to £ cov does not depend on the 
signature, and can therefore be done cither before or after 
the Wick rotation. 
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Canonical Approach. Given a particular folation of 
a globally hyperbolic Lorentzian manifold (with bound- 
ary) into Cauchy surfaces 11(f), we will assume that there 
is some way to define a canonical theory of free electro- 
magnetic fields, in terms of the on-shell degrees of free- 
dom on each time £(t), modulo gauge transformations. 
Rather than constructing this theory, we will take a more 
axiomatic approach and list the properties which we ex- 
pect it to have. 

In the canonical approach, the evolution operator 
£' can (Ei — > £2) must be a manifestly unitary evolution 
operator. (On a non-static background, the evolution op- 
erator £ is actually a map between two different Hilbert 
spaces. When we call such an evolution operator "uni- 
tary" , we really mean that it is an isomorphism between 
these two Hilbert spaces.) However, we do not assume 
that £ can = £ cov , or even that £ can is necessarily covariant 
under a change of foliation. 

Instead, we assume that the canonical and covariant 
theories both satisfy the classical linear field equations, 
interpreted as operator equations. Using these equations 
of motion, we can fix the evolution of the wavefunction 
from Si to £2 U P to an overall numerical factor X which 
can depend only on the geometry of M: 

(M, A aM ) = Z can {M, Aom)X(M). (49) 

B. Argument for unitarity 

We will now argue that £ cov is unitary in Lorentz sig- 
nature. The normal method for proving unitarity is to 
show that the covariant theory is equivalent to the canon- 
ical formulation using a reduced phase space. However, 
according to Ref. Q these two formulations arc not in 
general equivalent. Nevertheless, it is still extremely use- 
ful to know that a manifestly unitary canonical formula- 
tion of Maxwell exists, because we can use Eq. (|49|) to 
rewrite the covariant evolution operator as 

£ cov = X£ can = \X\U, (50) 

where X is a complex number, and U is some unitary 
evolution rule (i.e. an isometry). In the second equality 
we have performed a polar decomposition of X, and ab- 
sorbed the phase into U. After properly renormalizing 
ultraviolet divergences, both evolution operators should 
be finite, so X / 0. Hence one can define a unitary 
evolution operator as: 



Although we have used £ can to prove the existence of U, 
once we know it exists it can be defined using £ cov alone, 
by choosing the one and only positive real number \X\ 
that makes the right-hand-side unitary. So by Curie's 
principle (unique solutions to symmetric problems have 
symmetric solutions), U and |X| are themselves covariant 
and T- invariant. 



We will now show that |A| = 1, i.e. £ cov is unitary. 
First we observe that |A| by construction depends only 
on the geometry, not the state. Next we exploit the glu- 
ing property of Eq. (TTJ): 

£ cov (M! U M a ) = £ cov (M 1 )£ cov (M 2 ). (52) 

Because the product t/(Mi)[/(M 2 ) of two unitary oper- 
ators is unitary, and U{M\ U M 2 ) is the unique unitary 
operator equal to a real number times E cov (Mi U M 2 ), it 
follows that 

U(Mi U M 2 ) = U(M 1 )U(M 2 ). (53) 

Hence \X\ is also local: 

\X\(Mi U M a ) - |A|(Mi)|A|(M 2 ). (54) 

Using these relations, we can chop up the manifold into 
tiny pieces, in order to show that \X\ is given by a local 
action: 

\X\ = e iS M. (55) 

Here the action S is imaginary, and depends only on the 
metric g a b- This action might include boundary terms at 
dAI = Ei + £2- Temporarily ignoring these boundary 
terms, we can expand out the action in derivatives of the 
metric to obtain 

S = i d D x^g[ a + bR+ (cR 2 + dR 2 ab + eR 2 abcd ) + ...]. 

J M 

(56) 

When this action is Wick rotated back to the Euclidean 
signature, it remains purely imaginary. This is because 
each power of the Riemann tensor has an even number 
of time-derivatives, while the i coming from rotating the 
volume form \J—g = iJg turns e lS into e~ s . It follows 
that such terms cannot exist by T-invariance. (Although 
this fact could be deduced from the Lorentz signature 
theory, it is easiest to diagnose violations of T in Eu- 
clidean signature since there T consists of complex con- 
jugation only.) Therefore the covariant partition function 
is unitary up to a local boundary action. 

Terms in the boundary action could depend on the 
metric on the boundary, as well as its normal deriva- 
tives. However, it turns out that no boundary term can 
be found having the necessary properties. We start by 
mapping each of the initial and final Hilbert spaces into a 
specific Hilbert space H described by an orthonormal ba- 
sis of states (tpi, -02, "03--O- This allows U to be regarded 
as a unitary operator in H. In this basis, the inner prod- 
uct ($, \P) can be written as the identity matrix. Since 
raising and lowering with respect to the identity matrix 
is trivial, the boundary terms must take the same form 
for the "initial" slice £1 and the "final" slice £2. 

The absence of allowed boundary terms is easiest to 
see if lnlX I is rotated back to Euclidean signature. In 
order to satisfy the gluing property, the boundary term 
would have to contain an odd number of normal deriva- 
tives. (Terms such as Vol(dM) containing an even num- 
ber of derivatives would not cancel when gluing together 
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two regions.) T-reversal symmetry says the boundary 
term has to be real. Now an odd, real boundary term in 
Euclidean signature Wick rotates into a contribution to 
the Lorentzian ln|X| which is odd and imaginary. But 
ln|X| = iS is real in Lorentz signature by construction 
(and in any case, an imaginary contribution to iS would 
not lead to a violation of unitarity.) 

Since no contributions to \X\ are possible, it follows 
that £ cov is unitary. This is true even if it disagrees with 
£can on nonstatic manifolds. 



VI. DISCUSSION 

In section mi we defined the covariant electromagnetic 
partition function on spatially compact manifolds, taking 
into account issues involving zero modes of the field and 
gauge symmetry, as well as nontrivial U(l) bundles and 
large U(l) gauge transformations. 

Section IIIII showed how to implement the Faddeev- 
Popov trick on compact Euclidean manifolds. In Eq. (|l?p 
we found several kinds of departures from the usual "vec- 
tor minus two scalars" partition function. In addition to 
summing over nontrivial bundles, modding out by large 
gauge transfromations (i.e. Gribov copies), and omitting 
ghost zero modes, there is also a nontrivial normalization 
factor (q / [i^) y* 1 £ / 2irV , depending on the charge q, the 
t'Hooft gauge parameter £, the spacetime volume V, and 
the dimensionful parameter fi a which defines the measure 
on the space of gauge-transformations. 

Because the normalization factor depends on the ge- 
ometry through V, it contributes to the vacuum stress- 
energy tensor, and to the entropy of e.g. de Sitter 
space. In electromagnetism with one spacetime dimen- 
sion (which has no degrees of freedom at all) , this factor 
is necessary to obtain a trivial partition function. 

We also showed that the covariant approach is local 
and unitary, although the Faddcev-Popov trick obscures 
these properties. In our proof of unitarity ( section fV*)) . 
we relied on the existence of a canonical formulation of 
electromagnetism, but we did not explicitly construct the 
canonical formulation, nor did we prove that the two ap- 
proaches are generally equivalent. 

We were able to explicitly compare to the canonical 
formulation for two special categories of Euclidean man- 
ifolds (section IIV|) . One was D = 2 electromagnetism 
(where there are no local degrees of freedom), treated 
in [HI • The second was ultrastatic manifolds of the 
form S'xS, where the canonical partition function takes 
the form tr e~@ H summed over physical modes. In both 
of these cases, we found an agreement with the covariant 
approach. Although we have not done the calculation, 
we expect that the two approaches will always agree for 
static, horizonless manifolds. 

However, we expect that serious complications will 
arise for D > 2 time-dependent manifolds. In these 
cases, the two ghosts do not in general cancel out the 
temporal and longitudinal modes of [H, . In other 



words, the unphysical modes give rise to physical effects, 
which arc hard to explain in terms of the canonical de- 
grees of freedom. Similarly, in the case of manifolds with 
horizons, there are effects arising from the fact that vec- 
tors and scalars have different boundary conditions at 
the horizon. This seems to be related to the mysterious 
"contact term" discovered by Kabat [lH , who calculated 
that Maxwell fields contribute a negative divergence to 
the horizon entropy. This contact term cannot come from 
the entanglement entropy of the canonical modes, which 
is intrinsically positive. A similar discrepancy appears 
for gravitons [lj, HU • 

In fact, the original motivation of this work was to in- 
vestigate the validity of this contact term on compact 
manifolds (where one need not worry about complica- 
tions due to a boundary at infinity). In an earlier article 
[Toj , we showed that the contact term is absent in D = 2 
for canonically quantized electromagnetism. This article 
confirms that the covariant D — 2 theory is equivalent, 
so the contact term is absent there as well. The case of 
higher dimensions will be treated in a future article. 

Nevertheless, quite aside from the contact term, it 
would be disturbing if the covariant and canonical for- 
mulations gave different answers to physical questions. 
This article goes halfway towards resolving the issue, by 
showing how to properly treat the covariant formulation. 
Further progress on the discrepancy will require a better 
understanding of the canonical approach. 

It is not surprising that little work has been done on 
canonical Maxwell fields on nonstatic spacetimcs, be- 
cause the required conceptual apparatus is technically 
difficult. For each time slice, one has to construct a 
Hilbcrt space based on physical modes only. In order to 
write down the dynamics in terms of a time-dependent 
Hamiltonian, one must first select an isomorphism be- 
tween the Hilbert spaces at each moment of time (on a 
static manifold one can use time translation symmetry 
to do this, but on a nonstatic manifold one has to make 
some arbitrary choices). At the end, one has to prove 
covariance, i.e. that the dynamics would have been the 
same regardless of the choice of slicing. 

In Vassilevich's approach to canonical QFT, one slices 
up the spacetime, identifies the "physical modes" which 
would contribute to the canonical analysis, and then 
quantizes these modes using covariant methods Q . Vas- 
silevich found a noncovariant result, which disagrees with 
the covariant approach for e.g. Euclidean de Sitter space 
S 4 or S 2 . However, Ref. [|[ did not include the above 
normalization factor in his canonical approach. It also 
used the wrong measure in the path integral. (In the 
D = 2 case, the path integral measure on an expanding 
spacetime is a function of time. The correct measure is 
easiest to calculate if one uses the integrated s pac etime 
volume as the "time" parameter, as in Ref. llOj.) It 
may be that a more sophisticated version of the canoni- 
cal approach might give an answer which agrees with the 
covariant approach. 

The analysis of section IVl shows that the canonical and 
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covariant approaches agree up to a geometry-dependent 
phase factor. Therefore, so long as one is willing to 
insert such phases as a "fudge factor", the canonical 
method can be made to work. The question remains 
whether these phases can be derived elegantly in a canon- 
ical framework. 

In Wald's approach to canonical QFT in curved space- 
time [lH, one identifies the physical modes of the field 
in a slice-independent way, and quantizes this using the 
symplcctic structure. This approach is manifestly co- 
variant, but Wald was not able to directly calculate the 
stress-energy tensor. Instead he lists some axioms (in- 
cluding covariance) which the stress-energy tensor should 
satisfy, and used this to pin down their correct form. 

Assuming this approach can be applied to gauge fields 
(which Ref. [l6| claims, but does not work out in de- 
tail), it should be possible to derive the evolution op- 



erator £ can , use this to prove that £ cov is unitary, and 
then show that the covariant stress-energy tensor obeys 
Wald's axioms. Unfortunately, this roundabout method 
obscures the physical origin of the effects. It would be 
good to have a more direct way to derive the stress-energy 
tensor from the canonical formulation. If our treatment 
of the covariant calculation is any indication, this will 
probably require being very careful when handling the 
gauge symmetries. 
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